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Abstract. We present an analysis of the stellar mass growth over the last 10 Gyrs (z < 2) using a unique large 
sample of galaxies selected at 3.6/im. We have assembled accurate photometric and spectroscopic redshifts for ~ 
21200 and 1500 galaxies, respectively, with F(3.6pim) > 9.0/iJy by combining data from Spitzer-SWIRE IRAC, 
the VIMOS VLT Deep Survey (VVDS), UKIDSS and very deep optical CFHTLS photometry. We split our sample 
into quiescent (red) and active (blue) galaxies on the basis of an SED fitting procedure that we have compared with 
the strong rest-frame color bimodality (NUV — t')abs- The present sample contains ~ 4400 quiescent galaxies. 
Our measurements of the K-rest frame luminosity function and luminosity density evolution support the idea that 
a large fraction of galaxies is already assembled at z ~ 1.2, with almost 80% and 50% of the active and quiescent 
populations already in place, respectively. 

Based on the analysis of the evolution of the stellar mass-to-light ratio (in K-band) for the spectroscopic sub- 
sample, we derive the stellar mass density for the entire sample. We find that the global evolution of the stellar mass 
density is well reproduced by the star formation rate derived from UV based measurements when an appropriate 
dust correction is applied, which supports the idea of an initial mass function that is on average universal. 
Over the last 8Gyrs (z < 1.2) we observe that the stellar mass density of the active population shows a modest 
mass growth rate (p ~ 0.005(±0.005) M Q /Mpc 3 /j/r), consistent with a constant stellar mass density, p^ cUve ~ 
3.1 10 8 M Q /Mpc 3 . In contrast, an increase by a factor of ~ 2 for the quiescent population over the same timescale 
is observed. As a consequence, the growth of the stellar mass in the quiescent population must be due to the shutoff 
of star formation in active galaxies that migrate into the quiescent population. We estimate this stellar mass flux 
to be pa^q ~ 0.017(±0.004)Mq /Mpc /yr, which balances the major fraction of new stars born according to our 
best SFR estimate (p = 0.025(±0.003)M Q /Mpc 3 /yr). 

From z — 2 to z = 1.2, we observe a major build-up of the quiescent population with an increase by a factor 
of ~10 in stellar mass (a mass growth rate of ~ O.O63M0 /Mpc 3 /j/r). This rapid evolution suggests that we are 
observing the epoch when, for the first time in the history of the universe, an increasing fraction of galaxies end 
their star formation activity and start to build up the red sequence. 
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1. Introduction 

The strong decline of the cosmic star formation rate (SFR) 
since z ~ 1 is now well established (Schiminovich et al., 
2005; Hopkins et al., 2006), and has been shown to be ac- 
companied by a decrease of faint star forming galaxies and 
the decline of luminous ultra-violet galaxies (Lilly et al. , 
1996; Arnouts et al., 2005). One fundamental question is 
to understand the link between the global decline of the 
star formation rate and the evolution of the mass assem- 
bly. Inspection of the specific star formation rate (SFR per 
unit of stellar mass) in individual galaxies reveals that the 
preferred site of star formation activity has migrated from 
massive systems at high z to low mass systems at low-z. 
This is usually referred to as the downsizing effect (Cowie 
et al., 1996; Juneau et al., 2005). Evidences of such an 
effect have been seen in the fundamental plane relation 
of elliptical galaxies (Treu et al., 2005) as well as in the 
analysis of spectra of local galaxies (Heavens et al., 2004; 
Kauffmann et al., 2003) where both analyses show that 
massive galaxies have formed their stars earlier than less 
massive ones. Clustering properties of star-forming galax- 
ies at high and low redshifts reveal that the bulk of star 
formation activity has migrated from massive dark mat- 
ter halos (DMH) at high z to low mass DMH at low z, 
generalizing the downsizing effect from stellar mass to the 
dark matter (Heinis et al., 2007). 

While the ACDM hierarchical scenario successfully de- 
scribes the clustering properties of galaxies on large scales 
(Mo and White, 1996 ; Springcl et al., 2006), the quench- 
ing of star formation in massive systems is less well un- 
derstood. It relies on the complex physics of baryons and 
there is not yet enough constraints on the mechanisms 
involved in the regulation of star formation. These mech- 
anisms may act on galactic scales like AGN, Supcrnovae 
feedback (Benson et al., 2003, Croton et al., 2006) or on 
large scales via merging or gas heating in dense environ- 
ments (Naganimc et al., 2001; Yoshikawa et al., 2001). 
For instance, to reproduce the properties of local mas- 
sive elliptical galaxies, de Lucia et al. (2006) show that 
the stars could form first at high redshift in sub-galactic 
units, while the galaxies would continue their stellar mass 
assembly through merging onto low redshift. 
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In this context, the most stringent constraint for 
the semi-analytical models is the direct observation of 
the number density of massive galaxies at high redshift. 
Evidences of the existence of massive galaxies at high red- 
shift are becoming numerous. An important population of 
EROs at z > 1 has been discovered (K20 survey: Cimatti 
et al., 2002), with roughly half of them being old stel- 
lar systems (Cimatti et al., 2004). They could be the al- 
ready assembled progenitors of local massive ellipticals, 
as supported by their clustering signals and number den- 
sity (Daddi et al., 2002; Arnouts, 2003). The measures of 
the evolution of the stellar mass density from NIR surveys 
found that half of the stellar mass is already in place at 
z < 1-1.5 (Fontana et al., 2003; 2004; Pozzetti et al., 
2003; Drory et al., 2004; 2005; Caputi et al., 2005) mov- 
ing the formation epoch close to the peak of the cosmic 
SFR. The analysis of the galaxy mass function (GMF) at 
different redshifts provides additional information on how 
the process of mass assembly acts on different mass scales. 
Such analysis is becoming common among the recent sur- 
veys (Fontana et al., 2004, 2006; Drory et al., 2004, 2005; 
Bundy et al., 2005; Borch et al., 2006; Franceschini et al., 
2006, Pozzetti et al., 2007). While the different studies do 
not necessarily agree with each other, it emerges that the 
most massive galaxies have undergone less evolution than 
less massive ones since z ~ 1.0 and 50 to 70% are already 
in place at such a redshift, revealing a faster assembly for 
the most massive systems. 

An alternative approach to measure the density of 
massive galaxies at high redshift is to investigate the evo- 
lution of the local elliptical galaxies that dominate the 
massive end of the GMF (Bell et al., 2003; Baldry et al., 
2004, Bundy et al., 2006, Cirasuolo et al., 2006). This has 
been done by following the redshift evolution of the red 
sequence, defined by the (U — V) or (U — B) colors, from 
today up to z ~ 1.1. Bell et al. (2004) and Faber et al. 
(2005) find that the density of red galaxies drops by a 
factor 4 or 2.5 respectively, suggesting that their stellar 
mass density roughly doubles, while at the same time the 
blue population shows little evolution. They interpret this 
evolution by the migration of blue galaxies that quenched 
their star formation and migrated into the red sequence. 
Bell et al. (2004) also reported a significant evolution of 
the most massive ellipticals. Because there is a shortage 
of massive blue galaxies able to produce those massive 
ellipticals, they introduce the idea of "dry" or purely stel- 
lar mergers (merging between red galaxies) as a possible 
scenario to produce these massive local ellipticals. 

Analysis from the VVDS survey by Zucca et al. (2006), 
who define the early-type class based on SED fitting, re- 
port a more modest decline (~40%) in the number density 
of this population up to z ~ 1.1, which is also consistent 
with the analysis by Ilbert et al. (2006) based on a mor- 
phologically selected elliptical sample in the CDF South. 
Similarly to VVDS results, Brown et al. (2007), who have 
selected ellipticals with same optical color criteria than 
Bell et al. (2004), in an area of 7 deg 2 , do not observe any 
evolution of comoving density (<&*) up to z — 1 and ac- 
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cording to the modest luminosity density increase (~36%), 
predicts that the stellar mass has roughly doubled since 
z ~ 1. Moreover, they observe a modest evolution of lumi- 
nous ellipticals (with L > AL*) with 80% of their stellar 
mass in place at z ~ 0.8, suggesting that "dry" mergers 
should not play a dominant role in the evolution of ellip- 
ticals over the last 8Gyrs. 

Complementary analysis by Bundy et al. (2006) find 
signatures of a downsizing effect with quenching affecting 
first the massive galaxies and then moving to lower mass 
systems. This can be interpreted as an anti-hierarchical 
process with massive early types ending first their mass 
assembly while the less massive ones are still assembling 
their mass (see also Cimatti et al., 2006). If such a process 
is confirmed, it would be a strong observational constraint 
for the models of galaxy formation within the standard 
paradigm. 

A large fraction of the surveys discussed above are lim- 
ited to z < 1 — 1.2, due to optical selection, while the 
deepest surveys are still covering modest sizes in partic- 
ular when the sample is splitted by galaxy types. The 
impact of cosmic variance on results coming from small 
fields remains a major problem. Even in the three fields 
covered by COMB017, it still appears to play a signifi- 
cant role (Bell et al., 2004; Somerville et al., 2004). The 
selection of the samples varies from optical to IR band 
and the stellar mass is estimated from various approaches 
including SED fitting, rest-frame colors, with or without 
infrared information, introducing significant dispersions in 
the final estimates. Additionally, the selection of elliptical 
galaxies based on optical colors has been found to show a 
significant contamination by Sey2 AGN and star forming 
galaxies up to 40% (Franzetti et al., 2007). 

In this paper we explore the redshift evolution of the 
stellar mass assembly for a unique large sample selected 
in the observed Mid-Infrared, via the analysis of the K- 
rest luminosity functions, the luminosity density and the 
evolution of the mass to light ratio. We take advantage 
of a native IRAC selection that allows us to detect mas- 
sive red objects up to z ~ 2, and we isolate a large 
sample of red/quiescent galaxies according to their low 
level of star formation activity. We then discuss impli- 
cations for the stellar mass assembly of the active and 
quiescent populations. Throughout the paper we adopt 
the concordance cosmology il m = 0.3 and £l\ = 0.7 and 
Ho = 70km. s -1 Mpc" 1 . All magnitudes are in the AB sys- 
tem, and for convenience we adopt a Salpeter IMF unless 
otherwise specified. 

2. The data 

In this work we make use of the large amount of data col- 
lected in the VVDS-0226-04 field. The present sample is 
based on a 3.6/xm flux limited sample from the SWIRE 
survey which overlaps the deep multi-colour imaging sur- 
vey from the CFHTLS. The common area between the two 
surveys corresponds to 0.85 square degree. We complete 
this dataset by including the deep spectro-photometric 
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Fig. 1. Stars (triangles) and Galaxy (circles) number 
counts for the 3.6/im (red symbols) and A.bfim (green sym- 
bols) samples. Solid lines show the IRAC number counts 
from Fazio et al. (2005). 



data from the VVDS survey and the infrared data from 
the UKIDSS survey. 

2.1. The SWIRE sample 

The present sample corresponds to a 3.6/^m flux lim- 
ited sample with F(3.6/j,m)> 9/iJy (or toab(3.6) <21.5), 
based on the SWIRE survey (Lonsdale et al., 2003). The 
SWIRE photometry is based on the band-merged catalog 
including 3.6, 4.5, 5.6, 8.0/im and 24fim passband (Surace 
et al., 2005). with a typical 5cr depth of 5.0, 9.0, 43, 40 
and 311 /iJy respectively We use the flux measurements 
derived in 3 arcsec apertures for faint sources as suggested 
by Surace et al. (2005), while we adopt the adaptive aper- 
tures (Kron magnitude; Bertin and Arnouts, 1996) for the 
bright sources (m AB (3.6) < 19.5). In Figure [l] we show 
the IRAC 3.6 and 4.5/xm number counts. When comparing 
with the number counts from Fazio et al. (2005), the 3.6 
fim galaxy number counts are ^80% complete at 9/iJy. 
At this depth within the common area with CFHTLS, we 
have ~25500 sources. 

2.2. The CFHTLS data 

The deep multi-colour photometry (u*g'r'i'z') from the 
Canada-France-Hawaii Telescope Legacy Survey is based 
on the T0003 release (CFHTLS-D1). These data cover one 
square degree, with sub-arcsecond seeing in all bands and 
reach the limiting magnitudes (corresponding to 50% com- 
pleteness) of 26.6, 26.5, 26.0, 26.0 and 25.2 in u*, g', r', i', 
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z' respectively. A full description of the data will be pre- 
sented in a forthcoming paper by Mc Cracken et al. (in 
prep.). 

Thanks to the very deep optical data, almost all the IRAC 
sources have an optical counterpart, except a negligible 
fraction (~ 0.5%), due to unmatched positions within 
1.5arcsec. 

2.3. The VVDS data 

The VIMOS VLT Deep Survey consists of deep photome- 
try and spectroscopy (Le Fevre et al., 2005): 

• Deep B,V,R,I imaging with a depth (50% completeness) 
of 26.5, 26.2, 25.9, 25.0, respectively (Mc Cracken et al., 
2003). In addition, J and K observations with NTT-SOFI 
have been obtained at the depth (50% completeness) of 
24.2 and 23.8, respectively over 172 arcmin 2 (Iovino et 
al., 2005). 

• The first epoch VVDS spectroscopic sample is based 
on a randomly selected sample of ~ 9000 sources in the 
magnitude range 17 < I < 24.0 (Le Fevre et al., 2004). 
The spectroscopic area overlapping the SWIRE-CFHTLS 
data is 0.42 deg 2 and provides ~ 1500 secure spectra for 
sources with mAs(3.6/xm) < 21.5. 

2.4. The UKIDSS data 

Finally, we complete our dataset with the J and K pho- 
tometry from the UKIDSS Ultra Deep Survey (Lawrence 
et al, 2006) based on the DR1 release (Warren et al., 2007). 
These data reach a depth (5<r limits) of 22.5 and 22.0 in 
J and K, respectively. The overlap with SWIRE-CFHTLS 
is 0.55 and 0.76 deg 2 for J and K bands, respectively. We 
match the UDS sources with the optical data within a 1 
arcsec radius. Within the common areas, 89% (or 93%) of 
the 3.6 selected sources have a J (or K) flux measurement. 

2.5. The combined sample 

In Figure [2j we show the (i — 3.6/im) vs 3.6/im color- 
magnitude diagram for the whole sample (small dots). 
The galaxies and stars from the spectroscopic sample are 
shown as blue and green symbols respectively. We plot 
the cut-offs introduced in the spectroscopic sample due 
to the magnitude limit, Iab __ 24 (lower dashed line) 
and in the CFHTLS by adopting the 5<r detection limit, 
i' ~ 26 (upper dashed line) . The evolutionary tracks of an 
elliptical galaxy, formed at z = 6 and with an c-folding 
parameter r = O.lGyr (using the PEGASE model; Fioc 
et al., 1997) are shown as solid lines for different K band 
absolute magnitudes ranging from Kabs — ~ 25 to —21. 
While an optically selected sample with / < 24, can detect 
ellipticals brighter than Kabs —23.0 up to z ~ 1.2, a 
3.6^m selected sample with m(3.6/xm) < 21.5 can detect 
them up to z ~ 2 if very deep optical photometry is 
available. The present sample is therefore well adapted to 
investigate the evolution of old and/or massive galaxies 




18 19 30 Si 

3.6firn 

Fig. 2. Color-magnitude distribution with (i — 3.6//to) vs 
3.6/iTO. We show the spectroscopic sample for galaxies 
(filled blue symbols) and for stars (green stars symbols). 
The spectroscopic flux limit (/ < 24) and the 5cr detection 
limit for the CFHTLS (i' ~ 26) are shown as dashed lines 
(lower and upper line respectively). The solid lines show 
the behaviour of an elliptical galaxy (with r = O.lGyr 
and z/orm — 6) moving from z = 2.5 to z = (from 
red to blue color) and for different absolute K magnitudes 
between Kabs — — 25 (left solid line) and Kabs = —21 
(right solid line). The large circles denote the redshifts: 
z = 1.0, 2.0, while open squares are spaced by Sz = 0.2. 

up to z ~ 2, providing a major step with respect to 
previous analyses based on optical selection. 



3. Photometric redshifts 

3.1. The method 

We measure the photometric redshifts and we classify the 
whole population in galaxy/quasar/star based on the \ 2 
fitting analysis of the spectral energy distributions, using 
the photometric redshift code "Le Phare'Q In this work, 
we adopt a similar procedure to the one described by Ilbert 
et al. (2006) . We use empirical templates based on the four 
observed spectral types from Coleman et al. (1980) and 
add two starburst templates from Kinney et al. (1996). 
Templates are extrapolated into ultraviolet and infra-red 
wavelengths using the GISSEL synthetic models (Bruzual 
and Chariot, 2003) and we refine the set of SEDs with 
a linear interpolaton amongst the original SEDs. For Scd 
and later spectral types, we allow for different amounts 

1 http://serweb.oamp.fr/perso/arnouts/LE_PHARE.html 
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Fig. 3. Photometric vs Spectroscopic redshift comparison, 
for ~ 1400 secure VVDS spectroscopic redshifts. Open 
symbols refer to objects with internal photo-z uncertain- 
ties a{z phot ) > 0.3. 

of dust attenuation with a reddening excess E(B-V) vary- 
ing from to 0.6 and an interstellar extinction law from 
Prevot et al. (2004). The spectroscopic sample is used to 
perform a template optimization and estimate possible 
systematic shifts amongst the different passbands as intro- 
duced by Ilbert et al. (2006) . The SED fitting is performed 
from the u* to the 3.6 and 4.5/zm photometric passbands 
(the two latter passbands being still dominated by stellar 
light). All objects have at least 4 passbands with measured 
fluxes which ensures reliable SED fitting analysis. 
In Fig[3j we compare the photometric and spectroscopic 
redshifts for 1400 galaxies observed by the VVDS. We ob- 
tain an accuracy in the photo-z of cr[Az/(l + z)] ~ 0.031 
with no systematic shift and a small fraction of catas- 
trophic errors (1.5% with |Az| > 0.15(1 + z)). 
We consider as stars and QSOs (or AGN) the objects with 
a high SExtractor stcllarity index parameter in i' band 
(CLASS_STAR> 0.97) and for which either the star or 
the QSO template provide a best % 2 fitting. After remov- 
ing these stars and QSOs, we end up with a total sample 
of ^21200 galaxies out of 25500 sources. 



3.2. The Redshift distributions 

The redshift distributions are shown in Figure [4] for 
the whole sample with photometric redshifts (thick 
line), for the spectroscopic sub-sample (/ < 24, thin 
line) and for the photo-z 's of the spectroscopic objects 
(dashed line). The two latter distributions, both optically 
selected, are in excellent agreement, but do not show 



Fig. 4. The redshift distributions for the 3.6/im sample 
with F(3.6[im) > 9/iJy, for the full photo-z sample (solid 
thick line), the spectroscopic sub-sample (solid thin line) 
and the photo-z of spectroscopic objects (dashed line). All 
three curves are normalized to unity. 



the high redshift tail (up to z=2-2.5) observed for the 
3.6/im sample. We find that ~ 40% of galaxies with 
F(3.6/xm) > 9/zJy are at z > 1. A similar fraction, 
~ 35%, is obtained by Franceschini et al. (2006) for 
galaxies with F(3.6fxm) > 10/xJy in the GOODS-CDFS 
field. Note however that both samples rely partially on 
photometric redshifts. Since the high-z tail (z > 1.4) 
cannot be verified by our spectroscopic sample, we have 
investigated the reliability of the redshifts for this high-z 
population through color-color diagnostics. In Figure [5] 
we show the (g'-z') vs (z'-36/im) color-color diagram 
as an indicator for galaxies at z > 1.4. This is similar 
to the BzK diagnostic proposed by Daddi et al (2004). 
The locus of galaxies with z p h ot > 1.4 is well separated 
from the low redshift population as illustrated by the 
VVDS spectroscopic sample. This was also reported by 
Daddi et al. (2004) for the spectroscopic follow-up of BzK 
candidates. We conclude that our MIR selection at 3.6^m 
allows us to probe a significant population of galaxies in 
the redshift range 1 < z < 2 with no indication of any 
bias from the lack of spectroscopic redshifts. 

The spectroscopic redshift distribution also shows 
strong peaks at z ~ 0.35, 0.6 and 0.85. Two of them 
(z=0.35 and 0.85) are still observed in the photo-z dis- 
tribution although the area is two times larger and the 
distribution is somewhat smeared by the photo-z errors. 
These two peaks have been identified as large scale struc- 
tures extending through the entire field (Marinoni et al., 
2005; Adami et al., 2007) and we estimate the number den- 
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Fig. 5. Color-color diagram with (z' — 3.6/j,m)vs(g' — z') 
for 3.6/1 sample. The filled dots show the spectroscopic ob- 
jects while contours show the distributions of the photo- 
z sample for z < 1.4 (dashed-dotted green contours), 
z > 1.4 for the active/blue (long dashed blue contours) 
and stars (thin cyan contours). The quiescent/red popula- 
tions with Zphot > 1 -4 are shown as triangles and the open 
circles show dusty galaxies observed at 24/im (F(24^m) > 
300/iJy) and falling above the (NUV — t')abs criteria de- 
scribed in the text. The lines show a typical separation 
between the three classes in a similar way as the BzK 
selection. 

sity to be in excess by ~20-25% in these structures (based 
on a smoothed redshift distribution derived by applying 
a moving window with a width of Az — 0.2). Despite the 
large field of view, cosmic variance remains a significant 
source of uncertainty as previously pointed out by Bell et 
al. (2004) in the case of the COMBO 17 survey. Following 
the recipes described by Somerville et al. (2004), we esti- 
mate the cosmic variance from large scale structure fluc- 
tuations in the ACDM paradigm. In the low and high red- 
shift bins (with volumes of 2.7 and 12.5 10 5 Afpc 3 respec- 
tively), we expect the cosmic variance to produce density 
variations of 30% and 10% respectively, which appears to 
be in good agreement with our own estimates. We decide 
to account for this effect as an additional source of uncer- 
tainties in the measurements presented hereafter. 

3.3. Estimates of absolute magnitudes 

Throughout this paper we use absolute luminosity quan- 
tities. In practice, to derive the absolute magnitudes we 
use an adaptive filter method as introduced and tested on 
simulations by Ilbert et al. (2005, Appendix A). To re- 
duce the uncertainties in the k-corrections used to derive 



the absolute magnitude at X re st, we choose the apparent 
magnitude in the filter (X bs) that best matches the red- 
shifted Xrest- In the case where A 6 S = A rest * (1 + z), 
the k-correction is reduced to its redshift component 
(2.5 * log(l + z)) with no template dependency. 
All absolute magnitudes are derived according to this 
scheme except for the specific case of the K absolute mag- 
nitude (Kabs)- We always use use the 3.6/xm magnitudes 
for z > 0.4. This choice yields null k-correction at z ~ 0.6 
and requires small k-correction in highest redshifts with 
little template dependency. When available, at z < 0.4, 
we use the observed K band. 

In the next section we use the NUV and r' absolute lu- 
minosities. Over the redshift range 0.3 < z < 1.75 the 
rest-frame NUV (~ 2350A) and r' passbands correspond 
to observed wavelengths of 0.3/im < A < 0.65/im and 
0.8/ito < A < 1.7 fim respectively. We therefore need only 
small extrapolations for the NUV absolute magnitudes in 
the lowest z-bin. For the r' filter, the near infrared pho- 
tometry is required at high redshift, which is available for 
a large fraction of our sample. We note however, that even 
without near-infrared photometry the shape of the SED 
is always well constrained in the NIR domain, because all 
our objects have per definition a measurement at 3.6/xm. 
We therefore obtain a good estimate of the r' absolute 
magnitudes even at high redshift. 

We have derived the absolute magnitudes based on the 
empirical SEDs used to measure the photo-z. As a ex- 
ternal check, we have also measured the absolute mag- 
nitudes based on a library of SEDs from PEGASE2 code 
(Fioc et al., 1997). Thanks to the use of the adaptive filter 
method, the results from the two libraries are fully con- 
sistent with an r.m.s. dispersion for the (NUV — t'^abs 
color of <t[A(NUV — r')] ~ 0.2. This is small compared to 
the large dynamical range of this color (see Fig [6]) . 

3.4. The selection of quiescent galaxies 

We split our sample in two main classes in order to 
distinguish the red/quiescent and blue/active galaxies. 
To do so, we explore and compare two approaches, where 
one is based on the color bimodality and the other 
on an SED fitting procedure. For the color bimodality 
we consider the excellent separation provided by the 
(NUV — t')abs color (Figure[6] right panel) as suggested 
by recent GALEX results (Salim et al., 2005; Wyder et 
al., 2007). This is in contrast to the often used rest-frame 
(U - V)abs color suggested by Bell et al. (2004). 
One of the advantages of the (NUV — t')abs color 
is that the NUV passband is sensitive to a stellar 
population with a light-weighted age of < t >~ 10 8 yr 
while the r' passband is probing < t >> 10 9 yr (Martin 
et al., 2005), making the (NUV — r')ABS color an 
excellent indicator of the current over past star formation 
activity. Indeed, the GALEX-SDSS sample shows that 
(NUV — r')ABS is tightly correlated with the birthrate 
b parameter (b = SFR(t < 10 V)/ < SFR >) with 
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Fig. 6. Left panel: (NUV — t')abs color vs redshift for 
the quiescent/red (red contours) and active/blue (blue 
contours) populations. Color evolution for models with 
Zf = 3, 6 and r = 0.1 (dotted and solid lines resp.); Zf = 3 
and t = 1,1.5 (short and long dashed lines); r = 2, oo 
(dot-short dashed and dot-long dashed lines resp). Right 
panel: histogram for the (NUV — t')abs color (solid line). 
The dashed lines show the distributions of the two popu- 
lations as separated by SED templates (see text). 

(NUV — r')ABS > 4.3 being associated to galaxies with 
b < 0.1 and morphologically with de Vaucouleur profiles 
(Salim et al., 2005). 

To illustrate the redshift evolution in (NUV — r')ABS, 
we use the stellar synthesis model PEGASE2, with smooth 
star formation histories, described by different e-folding 
times and formation redshifts (r(Gyr), zj). In Fig [6] left 
panel, we show the evolution of (NUV-t')abs vs z for dif- 
ferent SF histories (r) and redshift of formation. Galaxies 
with t < 0.1 Gyr are in the red clump at all z, while mod- 
els with higher t move progressively into the red sequence 
at decreasing redshifts (z^0.7 and 0.4 for r = 1 and 1.5 
Gyr respectively). All models with r > 2Gyr have ongo- 
ing star formation that prevents them from reaching the 
red sequence unless their star formation is quenched pre- 
maturely. Therefore the (NUV — r') abs color bimodality 
(Fig [6j right panel) appears to be a good way to separate 
active and quiescent galaxies according to their birthrate 
parameter, with b < 0.1 for quiescent galaxies. 

In Figure [7J we show the galaxy distribution in the 
diagram (NUV — t')abs vs Kabs and the projected 
histograms for the different redshift bins. A strong bi- 
modality is observed at almost all redshifts providing a 
natural way to split the sample. We define an adaptive 




Fig. 7. (NUV — r')ABS color vs Kabs distribution (left 
side) as a function of redshift (as specified in each panel) . 
The two colored set of points correspond to the quies- 
cent and active galaxies as separated by SED templates. 
The solid lines denote the valley between the two popu- 
lations as seen in the projected histograms (right side). 
The green dotted histograms show the selection of red 
galaxies based on U — B rest- frame color (see text). The 
dashed histograms show the galaxies detected at 24^to 
(F(24^m) > 300/iJy) 

color cut-off located in the valley and characterized by: 
(NUV - r') AB s = -0.06 x (K ABS + 22) + b(z), where 
b(z) = 4.4 — 0.13 x z (solid lines). However, some uncer- 
tainties remain on the nature of the population in the red 
sequence due to the contamination by dusty star-forming 
galaxies. A reddening excess of E(B — V) = 0.2 produces 
a color reddening of A(NUV - r) A BS = 0.94 (or 1.04), 
assuming starburst (or SMC) extinction curve. We also 
show in Fig [7] the locations of those galaxies that are most 
likely dusty star- forming galaxies, as indicated by their de- 
tection at 24/im (with F(24[im) > 300/iJy). Between 5 to 
9% (depending on z) of this bright 24/im population lies 
above our adopted colour cut and is a source of contami- 
nation for our quiescent population. 

We then take a different approach, by classifying the 
galaxies according to the best-fit SED, on the basis of the 
six original templates (but allowing for dust extinction for 
the later types), in a similar way as Zucca et al. (2006). 
We consider as quiescent galaxies only those classified as 
elliptical and all other galaxies as active. 
In Figure [7] we compare the SED-based classification with 
the (NUV — r')ABS color distribution (quiescent as red 
points; active as blue points). To quantify the overlap be- 
tween the two different classifications, in each panel, we 
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Fig. 8. Examples of SEDs for a sample of high-z galax- 
ies with (NUV — t')abS above of the color criterion (see 
text). The fit is performed only in the wavelength domain 
0.3^m < A < 4.5/xm. The top four objects are best fitted 
by an old/passive elliptical, while the lower four galaxies 
are best fitted by dusty star-forming SEDs. 



measure the fraction of ellipticals in the sample above and 
below the color cut previously defined. We find that the 
SED-ellipticals dominate (with ~ 90%) the red peak pop- 
ulation up to z < 1.5 with only a small fraction falling 
below the cut. At z > 1.5, the fraction of ellipticals above 
the cut represents only 60% of the reddest objects while 
the remaining 40% are best classified as dusty star-forming 
galaxies. In Fig [8] we show examples of the two types of 
SEDs (Ell or dusty star-forming) for galaxies above the 
(NUV - r') cut and with z phot > 1.0. 
In Fig [5] we show the distribution of SED-elliptical galax- 
ies with z > 1.4 (red triangles). As for the BzK selection 
with a spectroscopic validation (Daddi et al., 2004), they 
are reasonably well separated from the active galaxies with 
z > 1.4 (blue contours). We also show the bright 24fim 
population with (NUV — t')abs m the red sequence and 
z > 1.4 (open circles). They are spread all over the upper 
part of the diagram. Amongst this population, ~50% of 
them has been correctly adjusted by dusty star-forming 
SEDs. This shows that SED fitting can get rid partially 
of the dusty galaxies falling in the red sequence. As a 
last test to quantify the residual contamination of the 
red/quiescent sample, we project our elliptical sample in 
the (J — K) vs (r' — K) color diagram. This is susceptible 
to distinguish the early-types from the dusty starbursts 
(Mannuzzi & Pozzetti, 2000; Daddi et al., 2004). We find 
that ~15% of the sample lies in the dusty starburst locus. 
If this color separation were to be efficient, this would 



mean that our quiescent sample may still be moderately 
contaminated by dusty galaxies. 

In conclusion, we have analyzed two methods to select 
the red/quiescent galaxies based either on the (NUV — 
t ')abs color or on the selection of quiescent galaxies ac- 
cording to SED fitting. We have shown that the SED- 
quiescent population represents the large majority of 
galaxies in the red peak of the bimodality distribution. 
We decide to adopt the SED-fitting classification because 
it allows to reduce the contamination of dusty starbursts 
in our red/quiescent sample. At z > 1.4, we estimate the 
residual contamination to be around ~15%. Only deeper 
Spitzer-MIPS observations could allow to better disentan- 
gle the two populations. However we note that the results 
discussed in the following sections are marginally affected 
by the chosen selection. 

We end up with a sample of 4425 quiescent and 16770 
active galaxies. 

4. The stellar mass to light ratio (M/L^) 

We now derive the evolution of the stellar mass to light 
ratio (M/Lk) with redshift for the active and quiescent 
populations. This will allow us to assess stellar mass quan- 
tities. Thanks to the large VVDS spectroscopic sample, we 
can measure how the M/Lk varies with redshift. In partic- 
ular, the use of spectroscopic redshifts allows us to avoid 
mass uncertainties relative to photo-z and we thus obtain 
accurate estimates of the scatter for these relations. 

We derive the stellar mass M* for each galaxy by fit- 
ting the multi-band photometry In the present work, we 
base ourselves on the Bruzual & Chariot (2003) stellar 
population model. We have used this model to construct 
a library of star formation histories that includes stochas- 
tic bursts (see Kauffmann et al. 2003 and Salim et al. 
2005 for similar libraries). Although this library uses a 
Chabrier (2003) IMF, we have decided to adopt a clas- 
sical Salpeter IMF throughout this paper. We therefore 
apply a systematic shift of +0.24 dex to the masses de- 
rived with a Chabrier IMF. In brief, our method uses the 
Baycsian probability distribution functions as described 
in Appendix A of Kauffmann et al. (2003) to derive pa- 
rameter estimates for each observed SED, including er- 
rors on each derived parameters, in this case the stellar 
mass. A full description of our mass estimation routines 
will appear in forthcoming papers (Lamareille et al. 2007, 
Walcher et al. 2007). For now we refer to previous works 
for discussion on the accuracy of the stellar mass esti- 
mates by SED fitting methods, in particular the choice of 
star formation histories with or without episodic bursts 
(Borch et al., 2006; Pozzetti et al, 2007) and the addition 
of Mid-IR bands (Fontana et al., 2006). We note, how- 
ever, that our stellar mass estimates may be affected by 
the established problems pertaining to the contribution 
of short-lived, luminous, infrared-bright stars in interme- 
diate age population (~ 1 Gyr) (Maraston et al. 2006). 
Our mass estimates would possibly need to be changed 
by 0.1 to 0.2 dex (see Pozzetti et al. 2007 for an in-depth 
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discussion). A final solution to this problem affecting all 
publicly available stellar population synthesis models is 
however not imminent. 

fn Figure [9] (upper panel) , we show the redshift distri- 
bution of the K luminosity (Lk / Lqk) for the quiescent 
(large black circles) and active (small grey circles) popu- 
lations. The solid lines show the evolution of the charac- 



teristic luminosity, L* , derived below in Section 5.1 and 
reported in Table [l] The distribution of M/Lk as a func- 
tion of redshift for the quiescent and active samples are 
shown in the middle and lower panels respectively. Before 
deriving the M/Lk versus redshift relations from the spec- 
troscopic sample we have tested several subsets: The first 
one includes only objects brighter than Lk > 10 11 Lqk 
to define an unbiased luminosity sample up to z = 1.5. 
The second one considers objects with a luminosity be- 
tween 0.4L*(z) < Lk(z) < 2.5L*(z), allowing to probe 
the evolution of the M/Lk for galaxies around L*(z) which 
are the main contributors to the total luminosity density 
with 70%, 45% and 55% for the quiescent, active and total 
samples, respectively. Because the spectroscopic sample is 
optically selected (/ < 24) , we have also considered a sub- 
sample with m3.6 < 20.5 (which means one magnitude 
brighter than the full sample) to insure a fair sampling of 
the {%'— 3.6/j.m) colors at all redshifts. All these subsamples 
are found to provide consistent results within the param- 
eter's uncertainties, therefore we only report the values 
for the three samples with the selection centered around 
their respective L*(z). For the three samples, we fit the 
redshift evolution of the M/L ratio by a power-law defined 
as log < M/Lk >= ax Z + b, with a global rms given by 
a. 

For the global sample, based on 999 galaxies: 

a = -0.30 ± 0.03, b = 0.03 ± 0.03, a = 0.22 (1) 

For the active/blue sample, based on 753 galaxies: 

a = -0.27 ± 0.03, b = -0.05 ± 0.03, a = 0.21 (2) 

For the quiescent/red sample , based on 298 galaxies: 



= -0.18 ± 0.04, b = +0.07 ± 0.04, a = 0.15 



(3) 



The fits to the quiescent and active samples are shown as 
solid lines in Fig [9j All subsamples show a similar trend 
with a gradual decline of the < M/L K > by a factor ~ 1.5 
to 2 up to z = 1. Similar results were obtained by Drory 
et al. (2004) for samples selected with different mass lim- 
its, based on the NIR MUNICS survey, and are shown 
in Fig [9] (long dashed lines) The behavior of the M/L is 
mainly driven by the stellar activity. The youngest and 
bluest stellar population have the smallest M/L \ (Bell and 
de Jong, 2001; Drory et al., 2004). Similarly, the decline 
of the M/L\ with redshift reflects the well established de- 
crease of the specific star formation rate (downsizing ef- 
fect; Cowie et al., 1996: for a galaxy with same mass, the 
SF activity is higher in the past). The smaller decline of 
M /Lk with redshift observed for the quiescent sample re- 
flects an older mean age of the stellar population, pushing 
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Fig. 9. Upper panel: distribution of the K luminosity 
{Lk I Lqk) for the spectroscopic sample (quiescent galax- 
ies: black filled circles; active galaxies: grey small circles). 
The characteristic luminosity L* for the Red/quiescent 
and Active/blue samples are shown as solid lines. Lower 
panels: Behavior of the Mass to light ratio (M/Lk) for the 
quiescent (middle panel) and active (lower panel) samples 
with the best fit shown as solid lines. Comparison with re- 
sults from Drory et al. (2004; long dashed lines) based 
on three different stellar mass cuts: M/Mq > 2.10 11 , 
10 11 < M/M© < 2.10 11 , and 4.10 10 < M/M & < 10 11 . 

the epoch of formation of the stellar component to high 
redshift. We find that the quiescent and active samples 
show a rather small scatter along this relation, with only 
0.15dex and 0.22dex scatter respectively. The main ori- 
gin for the scatter is the dependance of the M/Lk on the 
star formation history (SFH). The more complex SFH for 
the active population could be responsible for the larger 
scatter. 

5. The K band luminosity function and density 

5.1. The K rest- frame luminosity function 
5.1.1. K-LF measurements 

We measure the K rest-frame luminosity function (LF) by 
adopting the V max and STY estimators from the VVDS 
LF tool (ALF; Ilbert et 
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al. , 2004). In Figure 
show the LFs for the whole and spectroscopic samples 
and in Figure [ll] for the red/quiescent and blue/active 
samples. To measure the LFs, a weight has been applied 
as a function of apparent magnitude to account for the 
incompleteness in number counts. In the case of the spec- 
troscopic sample, we derive the LFs only for galaxies with 
m^B(3.6yu) < 20.5 and 0.2 < z < 1.2. This restriction is 
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Fig. 10. K-band rest-frame Luminosity functions, in 
different redshift slices, for the whole photo-z sample 
(Vmax- filled dark circles, STY: solid lines and shaded area 
based on slope uncertainties) and the spectroscopic sam- 
ple (V max : open circles). Also shown are the local LF from 
Kochanek et al. (2001; dotted line) and high-z LFs from 
Drory et al. (2003; green dahed-dot line), Caputi et al. 
(2005; blue long dashed line), Pozzetti et al. (2004; or- 
ange long dashed-dot line), Saracco et al. (2006; dark red 
long-short dahed line, a — —1.1), Cirasuolo et al. (2006; 
red dashed line). 

due to the I band selection effect (see Fig|2]). An additional 
weight is applied to the spectroscopic galaxies to take into 
account the sampling rate of the VVDS observations and 
its optical selection (as described in Ilbert et al., 2005). 
In Table [T] we report the values of the STY parameters 
for the different samples. The errorbars on the parame- 
ters refer to the Poisson errors and slope uncertainties. 
For the global and blue samples, we have fixed the faint- 
end slopes of the STY estimator at low and high red- 
shifts to the values observed between 0.6 < z < 1.0, 
where we have the best constraints in the bright and faint- 
end simultaneously. This choice is consistent with local 
(Kochanek et al. , 2001) and high redshift analysis (Drory 
et al., 2004; Caputi et al., 2005; Cirasuolo et al., 2006). 
For the red sample, we adopt a variable slope, although 
the low and high-z slopes are less well constrained. At 
low-z (z < 0.4) , we use a = —0.6, consistent with the 
slope observed by Bell et al. (2004) from a local ellip- 
tical sample. At z > 1, we adopt a gradual flattening 
of the slope that best reflects the evolution of the V max 
estimator. However the depth of the current data does 
not allow any statement on the reality of this apparent 
flattening and deeper observations are required to con- 
firm it. We note that the choice of the slope does not 



affect the discussion in next sections because changes to 
the estimates of the luminosity densities are marginal 
(p L (a) - p L (a = -0.3) = -0.01, -0.02, -0.04dex for 
< z >=1.1, 1.35, 1.75, resp.). 

5.1.2. Additional source of uncertainties 

Additional sources of uncertainties on the LFs come from 
the photo-z and cosmic variance. We have estimated the 
uncertainties relative to the photo-zs from 100 mock sam- 
ples based on the redshift probability distribution (PdZ) 
of each source (as in Bolzonella et al., 2000) and by re- 
computing the type separation and luminosity functions. 
Considering only the variations of <&*, we find uncertain- 
ties in the density to be between 5% to 12% for the total 
and active samples, depending on the redshift bin. For the 
quiescent sample, the uncertainties vary between 7% and 
28%. This method allows to account simultaneously for 
the large photo-z uncertainties of the ellipticals at high- 
z and to estimate the stability of the quiescent vs active 
separation at high redshift. We report these uncertainties 
in Table [l]((d$/$) Pdz)- Except for the quiescent sample 
at high z, the photo-z uncertainties are not a dominant 
source of errors in the present work. 

It is more difficult to evaluate the effect of cosmic vari- 
ance. As a first attempt we have measured the global 
LFs in 10 non overlapping fields covering the whole area 
and find variations of <&* of between 0.07dex to 0.2dex. 
However, the fact that the 10 fields are located in a sin- 
gle area of the sky prevents us from using them as cos- 
mic variance uncertainties over the whole survey (adopt- 
ing a simple rescaling as %/l0). Therefore, as mentioned 
in Section [2j we adopt the formal approach discussed by 
Somervillc et al. (2004). The cosmic variance is estimated 
as a = baoM, where (Jdm is the variance of the dark mat- 
ter that depends on the comoving volume (their Figure 3, 
right panel) and b is the bias relative to each specific sam- 
ple that depends on the comoving number density (their 
Fig 3, left panel). Theses uncertainties are reported in 
Table [I] ((d$/$)cv)- Cosmic variance remains the domi- 
nant source of uncertainty for the whole and active sam- 
ples at all redshift, while for the quiescent population the 
photo-z uncertainty is the major source of uncertainty at 
z > 1.35. 

5.1.3. Evolution of the K-band LFs and comparison 
with other studies 

In Fig [lOj we compare the global LFs with other stud- 
ies. Beyond the fluctuations of <I>* due to the presence 
of large scale structures, our LFs appear slightly brighter 
than previous studies, by about 10%. This is however con- 
sistent with the expected calibration uncertainty of the 
SWIRE-IRAC flux (Surace et al., 2005). Comparing K* 
to the local measurement by Kochanek et al. (2001), we 
observe a brightening of AK* ~ — 0.8( — 1.0) ± 0.2 from 
z — to z — 1.5(2). Regarding the comoving density pa- 
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Fig. 11. K-rcst luminosity functions for the blue/active 
(blue squares and solid lines) and the red/quiescent (red 
triangles and solid lines) samples. The thin black lines re- 
fer to the global sample from Fig [lOj We compare our LFs 
with the red and blue samples from the UKIDSS survey 
(Cirasuolo et al., 2006; red and blue dashed lines). 



rameter (<&*) we observe a decline by a factor 1.7±0.2 up 
to z = 1.2 and 2.5±0.5 up to z — 2.0 as compared to local 
estimate. This mild luminosity brightening and modest 
density decline agree with previous studies (Caputi et al., 
2005; Saracco et al., 2006, Cirasuolo et al., 2006). 



In Fig 11 we compare the LF per types with our global 
LF and results from UKIDSS based on (U — B) rest color 
selection (Cirasuolo et al. , 2006) . Our active sample shows 
a similar brightening than the global sample. The density, 
appears stable between 0.2 < z < 1.5 and starts to 
decline at higher z. For the quiescent sample a different 
behaviour is observed. Up to z — 1.2, the density declines 
by a factor 21J 7 , followed by a sharp drop, by a factor 
16tg, between 1.2 < z < 2. While the behaviour of the 
active and quiescent samples agrees with UKIDSS results 
on blue and red galaxies (concerning faint end slopes and 
global evolution), we obtain a different ratio between the 
blue/active and red/quiescent populations and they do 
not observe the strong decline for the red population at 
z > 1.2. We suggest two differences in the selection crite- 
rion for red/quiescent galaxies as possible source to this 
discrepancy: 

— In Fig [7] we show the (NUV — r') color distribution 
(green dotted histograms) for galaxies selected on the 
basis of the optical color (U — B), in a similar way 
as Cirasuolo et al.. While the Cirasuolo selection in- 
cludes our quiescent sample, it also includes galaxies 
with (NUV — r) as blue as 2.5 — 3. This means that the 



optical criterion yields a number of red galaxies that 
is larger by ^40 to 50% depending on redshift. 
— the modest decline of the red LF at z > 1.2 for 
UKIDSS can also takes its origin in the color crite- 
rion and in the inclusion of the larger fraction of red 
dusty starbursts at high-z. Regarding our selection of 
a quiescent sample, we were conservative in the sense 
that we required that the galaxies do not reveal signs of 
star formation activities according to their (NUV — r') 
color and by using multi-color SED fitting to exclude 
dusty starbursts. In our red population at 1.5 < z < 2, 
only 60% of galaxies is considered as quiescent, the rest 
being dusty starbursts. We expect this effect to be even 
stronger in the case of an optical color selection. 

The strong evolution between z = 2 and z = 1.2 of the 
LF parameters of the quiescent sample alone suggests that 
we are probing the epoch when an increasing number of 
galaxies stop their star formation activities and turn quies- 
cent. However, the LF parameters are strongly correlated 
with each others and their interpretation shaky. We there- 
fore now turn to discussing the evolution of the different 
samples via the measurements of integrated quantities (lu- 
minosity densities and stellar mass densities). 

5.2. The K band luminosity density 

A sensitive test for galaxy formation is the measurement 
of the luminosity density which provides the total amount 
of light emitted per unit of volume and is estimated as: 
PLk = f °° LQ(L)dL = $*L* K T(a + 2), where L is the 
gamma function and (<&*, L*, a) are the Schechter parame- 
ters listed in Table[l] The measurements of Pl k (expressed 
in solar unit with M^q = 5.14) are listed in Table [ll and 
plotted in Figure |12| The quoted errorbars include all the 
sources of errors discussed in previous sections. 
The total sample shows a luminosity density increase by 
a factor ~ 2 (starting from local values) up to z ~ 1 — 1.2 
followed by a fall off of a similar amount up to z = 2. This 
global behavior agrees, within the errorbars, with previous 
K band luminosity density studies and is well described by 
the fit proposed by Caputi et al. (2005, solid line). Looking 
at the SED selected samples, the active sample dominates 
the K band luminosity density at all z and follows a similar 
evolution as the global sample. For the quiescent sample, 
the luminosity density remains roughly flat up to z ~ 1.2 
then drops significantly. As discussed in previous section, 
this trend is less pronounced in the UKIDSS 's red sample 
measurements. 

To interpret the evolution of the luminosity density 
for the different samples, we need to set local references. 
Since no similar K-band based selection exists, we make 
use of existing optical catalogs with information that can 
be converted to yield a K band local reference. We have 
used the results reported by Driver et al. (2006) based 
on the Millennium Galaxy Catalog (MGC) where they 
provide the LFs in B band and stellar mass density (p*) 
for a large variety of classification, based on morphology, 
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Fig. 12. Evolution of the luminosity density, p^ K (ex- 
pressed in Lqk ), for the total (filled circles) , active 
(filled squares) and quiescent (filled triangles) populations. 
Results are compared with previous works as specified in 
the caption and with a pure luminosity evolution models 
(dashed lines). 

colors, SED fitting and spectral classes (77 parameter from 
2dF: Madgwick et al.; 2002). We adopt two criteria for 
the local reference that can reflect the definition of our 
quiescent sample: 

— As noted by Madgwick et al. (2003), the rj parameter 
is highly correlated with the birthrate b parameter. In 
particular their first class with rj < —1.4 corresponds 
to b < 0.1 and is thus similar to what we expect locally 
for our criterion. We adopt the rji class as representa- 
tive of our quiescent sample and define all other classes 
as belonging to the active sample as reported in Table 
2 of Driver et al. (2006) 

— The majority of galaxies in the (NUV — r') red se- 
quence observed by GALEX are dominated by de 
Vaucouleur profiles (Salim et al, 2005). We adopt the 
E/SO(red) morphological class of the MGC as re- 
ported in their Table 2 as another representation of 
our quiescent sample. 

Finally, we converted /?* to p^ K for each sample by 
adopting the Bell and de Jong (2001) relations between 
(B-R) and M/Lr-. The resulting ranges of local estimates 
are shown in the Fig[l2]as colored rectangles. We addition- 
ally plot the morphologically separated luminosity densi- 
ties from 2MASS (Kochanek et al. , 2001) as colored open 
stars. The low-z estimates agree reasonably well betwen 
themselves except for the quiescent sample from 2MASS. 
This is most likely due to the inclusion of the blue com- 
pact ellipticals in the sample of Kochanek et al. (2001). 



The MGC with its red E/S0 sample appears to be most 
suited to provide our local reference in this case. 

Using the local normalization, we can now compare 
the evolution of p^ K for pure luminosity evolution (PLE) 
models, i.e. with no merging involved. We normalize to 
the local density <S>*(z = 0) and use ALogpL K (z) = 
—OAAKabs(z)- We adopt the luminosity evolution as 
derived in Section |5.1| Specifically, for the total sam- 
ple we adopt the luminosity evolution from Caputi et 
al. (2005), for the active sample we adopt AKabs( z ) = 
— 0.65(±0.07) x Z, while for the quiescent sample we 
use two extrem PEGASE models that encircles the ob- 
served luminosity evolution and can be described by 
AK A bs(z) = -0.5| - 0.85 x Z. In Figll2j the PLE pre- 



dictions are shown with dashed lines and shaded areas. 
The evolution of the rest-frame K-band luminosity den- 
sity is consistent with the PLE model for the total and 
active samples up to z = 1.2, while the model overpre- 
dicts pl k at higher redshift. This shows that the number 
density of galaxies has to drop. For the quiescent sample, 
the PLE model fails to reproduce the observed trend at 
low and high redshifts. It predicts an increase in luminos- 
ity density by a factor 1.6 to 2.2, up to z ~ 1.2, while the 
observations suggest a modest increase by a factor rang- 
ing from 1.0 to 1.4. The disagreement with the PLE model 
at highest redshift is even more pronounced than for the 
active sample suggesting that the number density of the 
quiescent galaxies must drop even faster. 

6. The stellar mass density up to z=2 

We now derive the stellar mass density, p*, up to z=2. 
To that end we convert our K-band luminosity densi- 
ties (Table [T|), to stellar mass densities via the relation 
p*(z) = pL K (z)x < M/L^ > (z) and using the mass to 
light ratio equations determined in Section [4] 
Our measurements of p* for the three samples are shown 



in Figure 13 The errorbars account for Poisson, photo-z 
and cosmic variance uncertainties (as for p^) and an ad- 
ditional uncertainty of 0.05dex from M/L estimates (as- 
suming that the paramaters a and b in the M/L relation 
are un-correlated) . 

6.1. The evolution of the stellar mass density 

When moving back in time, the global population shows a 
small but regular decline up to z ~ 1.1 which accelerates 
at higher z. By comparing with local estimates, we find 
that the total stellar mass has decreased by roughly a 
factor of - 1.5, 2 and 4 up to z ~ 1.0,1.5,2.0. The 
compilation of previous surveys shows a large scatter 
by roughly a factor two and our estimates are located 
in the upper envelope. While the scatter can be in part 
due to mass estimates, for example the use or not of 
near IR data, the cosmic variance is most likely to be 
the dominant factor as discussed here and by Bell et al. 
(2003). 

We quantify the stellar mass evolution of the active 
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Fig. 13. Evolution of the stellar mass density as a function of cosmic time (assuming a Salpeter IMF). The total, active 
and quiescent stellar mass densities from this work are shown with large filled circles, blue squares and red triangles 
respectively. For reference, the right-hand axis gives the stellar density parameter. The integrated Star Formation Rates 
for different dust attenuation corrections (A-fuv = 1-1,1.3,1.7), based on the SFR derived by GALEX (Schiminovich 
et al., 2005), are plotted as solid lines, and the one from the compilation of Hopkins and Beacom as a dashed line. 
The dust corrected SFRs are shown in the inset. High z measurements from the literature for total samples have 
been splitted between analysis based on optical information only (grey symbols) and including Near or Mid IR data 
(green symbols) for the mass estimates, with Optical : Brinchmann & Ellis (2000; pentagons); Cohen et al. (2002; 
losanges); Dickinson et al. (2003; squares); Gwyn et al. (2005; inclined triangles); Borch et al. (2006; crosses), with 
NIR data: Drory et al. (2004; open stars) ; Drory et al. (2005; up and down triangles); Fontana et al. (2006; losanges); 
Franceschini et al. (2006; ellipticals: red down triangles and global samples : open squares); Pozzetti et al. (2007; 
open circles); Abraham et al. (2007; ellipticals : red up triangles). Local values are from Kochanek et al. (2001) for 
morphologically selected ellipticals (red star) and spirals (blue star) and whole sample (black star); Cole et al. (2001; 
circle) and Driver et al. (2006; colored rectangles). For clarity we do not show errorbars for other surveys. 



and quiescent samples with a simple linear fit with 
redshift up to z — 1.2. Including local measurements, 
we get : ALogp^ cttve = -0.05(±0.09)z + 8.51(±0.04) 
and ALogp+ uiescent = -0.31(±0.07)z + 8.38(±0.02). The 
active population shows a modest evolution, consistent 
with no evolution with a mean value p^ ctlve = io 8 - 49±0 04 . 
This constancy of the active (blue) sequence has been 
pointed out by Borch et al. (2006) and by Martin et 



al. (2007) who derived a similar value. On the other 
hand, we observe that the stellar mass of the quiescent 
galaxies, pQ mescent ^ nas increased by a factor 2±0.3 
between z = 1.2 and z — 0. Between z = 2 and 
z — 1.2 the evolution in stellar mass is even stronger, 
it increases by a factor of ^10. This suggests two dif- 
ferent regimes in the build-up of the quiescent population. 
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6.2. Comparison to the integrated star formation rate 
estimates 

Thanks to the recent constraints on the cosmic star for- 
mation rate, we can compare the observed stellar mass 
density (p*) to the predictions for the cosmic star for- 
mation rate (p*) over time, since the two quantities are 
simply related as follows: 

p*(t) = (1 -R) f p*{t)dt (4) 
Jo 

where R is the recycling factor (fraction of stellar mass 
released in the interstellar medium) which is assumed, in 
the case of a Salpeter IMF, to be R = 0.3 (Madau et al.; 
1998). For the SFR, p+(z), we use the uncorrected UV 
SFR as measured by Schiminovich et al. (2005) between 
< z < 3 that we have extended up to z ~ 6, based 
on measurements from Steidel et al. (1999), Bowens et 
al. (2006) and Bunker et al. (2004). To derive the dust 
corrected UV-SFR, we adopt a range of dust attenuation 
that varies between 1.1 < Apjjy < 1.7 which is consistent 
with local estimates (Buat et al., 2005) and with other 
dust-free SFR estimators up to z ~ 1 (Schiminovich et al., 
2005; Hopkins and Beacom, 2006). We also use the SFR 
estimate from the analysis of Hopkins and Beacom (2006) 
who compile most of the recent surveys from the UV to the 
Far-Infrared to derive a global, dust cor rected SFRQ The 
two SFR histories are shown in an inset in Fig [13] and are 
in good agreement within the range of dust attenuations 
we have considered. The result of integrating is shown 
as a shaded area. We emphasize that we did not impose 
beforehand that our integrated SFR be consistent with 
the stellar mass density at z = 0. Nevertheless, it provides 
a remarkable fit in the redshift range probed in this work 
and encloses well the local estimates (see also Rudnick et 
al., 2006). 

We find that our measurements are well represented by 
an attenuation value ranging from 1.1 < Apuv < 1-3) be- 
tween < z < 1.5. This constancy of the mean dust atten- 
uation up to z ~ 1.5, can be qualitatively interpreted as a 
trade-off between two competing factors: -1- star- forming 
galaxies at high-z have less dust attenuation than their 
local counterparts at the same total luminosity (Reddy 
et al., 2006, Xu et al., 2007) -2- star-forming galaxies at 
high-z are intrinsically more luminuous (Arnouts et al., 
2005; Le Floch et al., 2006). We can quantify these ef- 
fects by looking at the relations between the (Loust/Luv) 
ratio and the total luminosity (Ltot — Lp) USt + Lw), 
all quantities expressed as vh v . From the local relation 
(Martin et al. , 2005), a typical L* uv {z = 0) galaxy has a 
LDust/ Luv ~ 3.75, corresponding to Apuv ~ 1-3 (Buat 
et al., 2005). If we let this galaxy evolve according to the 
evolution of the luminosity functions (Lpuv ~ (1 + z) 2 5 ; 
Arnouts et al., 2005 and L FIR ~ (1 + z) 3 - 2 ; Le Floch 

2 We adopt their SFR based on the modified Salpeter IMF 
which is converted to a Salpeter IMF by applying +0.2dex in 
stellar mass. 



et al., 2006), at z ~ 1.35 we obtain a typical luminosity 
of log(LTOT I Lq) ~ 11.34. We can finally estimate the 
amount of dust attenuation for this galaxy based on the 
relation from Reddy et al. (2006, their Eq. 5) that holds in 
the redshift range 1 < z < 3. We obtain L Dust /L uv ~ 4.3, 
corresponding to Apuv ~ 1-4. This rough estimate shows 
that a typical L*{z) galaxy that contributes to the bulk 
of the SFR at all z shows, on average, a small evolution 
of the amount of dust attenuation, in agreement with the 
observation of the integrated dust corrected SFR in Fig[T3| 
To summarize, with simplified but not unrealistic es- 
timates of the dust correction, we observe that the two 
complementary quantities, the star formation rate and 
the stellar mass density, are remarkably consistent with 
each other. This result is encouraging, and may provide 
some more fundamental informations about the initial 
mass function (IMF). The observations of the two quan- 
tities explore two different parts of the IMF (dominated 
by massive stars, a few M©, for the SFR, and low mass 
stars, for stellar mass) and are extrapolated through the 
adopted IMF. The good agreement of these complemen- 
tary observations over the last 10 Gyrs, is a good support 
for an IMF that is on average universal, as already pointed 
out in previous work (e.g. Franceschini et al., 2006). We 
note that this remark is not specific to the Salpeter IMF 
adopted throughout this paper, but holds also for other 
IMFs suggested in the literature, as long as the quantities 
(p* , p+) are appropriately corrected for IMF effects. 

7. Discussion and Conclusions 

In this work, we have used a unique large sample of 
21200 galaxies selected at 3.6pm (based on the VVDS 
and SWIRE surveys) to investigate the evolution of the 
luminosity functions, luminosity densities, stellar mass to 
light relations and stellar mass densities up to z ~ 2. We 
have separated the active and quiescent galaxies based on 
an SED fitting procedure. We define as quiescent those 
galaxies which are best fit by an elliptical template. We 
have shown that this sample reproduces well the red se- 
quence of the color bimodality (NUV — r') (known to be 
a good separator between active and quiescent galaxies 
with b < 0.1) and at the same time minimizes the con- 
tamination by dusty starbursts. We find that the active 
and quiescent populations follow different behaviours. In 
particular, there is a clear transition between two regimes 
in the evolution of the quiescent population at a redshift 
of z = 1.2 (i.e. 8 Gyr ago). 

7.1. The Last 8 Gyrs 

The active population appears to be in place at z ~ 1.2 
with a small evolution over this time laps except the aging 
of its stellar population. In contrast the quiescent sequence 
shows a gradual increase with a doubling of its stellar 
mass in line with previous optical studies (Bell et al., 2003; 
Faber et al., 2006; Brown et al. (2006). 
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Table 1. K-LF parameters (a, M* , $*) and K luminosity densities (pl) for the whole, active and quiescent samples. 
Errors in M* , refer to Poisson errors and slope uncertainties. Additional uncertainties due to photometric redshift 
(PdZ) and cosmic variance (CV) are reported separately. Errors in pl include all of them. 



< z > 


# 


a 


M* 


$* 






PL 










(10~ 3 Mpc" 3 ) 






(lO 8 L 0if Mpc- 3 ) 


ALL 


0.3 


2180 


-1.1 ±0.2 


-22.84+'- 


4.19li9 


0.05 


0.23 


6.97 ± 2.2 


0.5 


2680 


-1.1 ±0.2 


-22.83 ±0.3 


3.50 ± 1.3 


0.05 


0.20 


5.77 ± 1.3 


0.7 


3336 


-1.1 ±0.2 


-22.96 ±0.2 


3.36 ± 1.1 


0.05 


0.16 


6.24 ± 1.3 


0.9 


4545 


-1.1 ±0.2 


-23.08 ±0.2 


4.22 ± 1.2 


0.05 


0.26 


8.75 ±2.5 


1.1 


3027 


-1.1 ±0.2 


-23.22 ±0.2 


2.75 ±0.7 


0.11 


0.13 


6.48 ± 1.4 


1.35 


3077 


-1.1 ±0.2 


-23.18 ±0.17 


2.56 ±0.5 


0.09 


0.12 


5.82 ± 1.5 


1.75 


2189 


-1.1 ±0.2 


-23.42 ±0.14 


1.47 ± 0.2 


0.07 


0.10 


4.17 ± 1.1 


active 


0.3 


1666 


-1.3 - ±0.2 


99 oq+U.43 
ii - OJ -0.S6 


2-38lij 


0.05 


0.27 


4.76 ± 1.7 


0.5 


2063 


-1.3 ±0.2 


-22.82 ±0.3 


2.22 ± 1.0 


0.05 


0.21 


4.40 ± 1.3 


0.7 


2569 


-1.3 ±0.2 


-22.95 ±0.2 


2.21 ± 0.9 


0.05 


0.17 


4.94 ± 1.3 


0.9 


3428 


-1.3 ±0.2 


-23.06 ±0.2 


2.82 ± 1.0 


0.05 


0.28 


6.97 ±2.7 


1.1 


2358 


-1.3 ±0.2 


-23.19 ±0.2 


1.99 ±0.6 


0.12 


0.14 


5.55 ± 1.8 


1.35 


2609 


-1.3 ±0.2 


-23.19 ±0.18 


2.08 ±0.5 


0.09 


0.13 


5.80 ±2.0 


1.75 


2081 


-1.3 ±0.2 


-23.49 ±0.15 


1.27 ± 0.3 


0.07 


0.11 


4.66 ± 1.8 


quiescent 


0.3 


514 


-0.6 ±0.2 


99 Q1 + u - a4 
^.a±_ _49 


l-78i u .f 


0.09 


0.33 


2.63 ± 1.1 


0.5 


617 


-0.3 ±0.2 


-22.55 ±0.19 


1.47 ±0.17 


0.07 


0.24 


1.60 ±0.4 


0.7 


767 


-0.3 ±0.2 


-22.73 ±0.15 


1.29 ±0.13 


0.09 


0.21 


1.65 ±0.4 


0.9 


1117 


-0.3 ±0.2 


-22.83 ±0.14 


1.58 ±0.14 


0.07 


0.34 


2.22 ±0.8 


1.1 


669 


0.0 ±0.3 


-22.86 ±0.17 


0.90 ±0.08 


0.12 


0.17 


1.42 ±0.3 


1.35 


468 


0.3 ±0.3 


-22.84 ± 0.16 


0.41 ±0.03 


0.21 


0.18 


0.74 ± 0.2 


1.75 


108 


0.6 ±0.3 


-23.27 ±0.15 


0.05 ±0.02 


0.28 


0.17 


0.17 ±0.05 



The lack of evolution of the stellar mass density for the 
active, p<^ ctlve j and the doubling for quiescent, p'i ule8cent j 
is surprising. The gradual increase of the global stellar 
mass density, over the last 8 Gyrs can be easily explained 
by the formation of new stars by star forming galaxies 
at a rate described by the cosmic star formation rate. 
Or, similarly, it is found that intermediate mass galax- 
ies (M/Mq ~ 10.5) have specific SFRs allowing them to 
increase their stellar mass by a factor two since z = 1 
(Juneau et al., 2005). On the other hand, the increase of 
^quiescent j g j n a pp aren ^ contradiction with the definition 
of quiescent galaxies that cannot produce more than 5% of 
additional stellar mass over the last 8Gyr (due to the ab- 
sence of star formation) . The most plausible explanation 

c ii r quiescent • r 

tor tne increase ot pi is a progressive migration ot 

galaxies from the active to the quiescent population at a 
similar rate than that with which new stars are formed. 
We quantify this evolution by a rough estimate of the stel- 
lar mass growth per unit of time defined as : 

PX - 7 7 _ , 7-T (5) 

^univ\ z l) L univ\ z h ) 

, where X refers to stellar mass of the considered sam- 
ple in the redshift range z ; < z < z h . We find that the 
stellar mass growth of the active population evolves very 
little with pAetive = 0.005 ± 0.003M e /Mpc 3 /yr, while 
the quiescent population has a much higher stellar mass 
growth with pQuieecent = 0.017 ± 0.004M o / Mpc 3 /yr. We 



can compare this value with the mass growth expected by 
integrating the star formation history over the same period 
of time. We find that p S FR = 0.025 ± 0.003M G /Mpc 3 /yr 
(assuming 1.1 < Apuv < 1-3). Under the assumption that 
the quiescent population has negligible residual star for- 
mation, its mass growth can be attributed to the mass 
flux of active galaxies moving into a quiescent mode 
(pA^Q = PQuiescent), and which appear to account for 
most of the global stellar mass growth derived from the 
SFR. 

Additional evidence for this scenario is also presented by 
Vergani et al. (2007) who split the VVDS spectroscopic 
sample in active and quiescent galaxies based on 4000 A 
break. In order to satisfy the constraints from the stellar 
mass functions split in active and quiescent galaxies, they 
find that the star formation in some massive blue galaxies 
must have been quenched, moving these galaxies into the 
"red sequence". 

7.2. The major epoch of build-up for quiescent 
galaxies: 

The present analysis allows us to extend previous work to 
the redshift range between 2 > z > 1.2 (3.2Gyr < T univ < 
bGyr). This appears to be a transition epoch, as evidenced 
by the strong increase in stellar mass, by a factor of ~10, 
of the quiescent population between 2 > z > 1.2, while 
the active population increases still by a factor ~2.5. This 



16 



Arnouts et al.: K-LFs and stellar mass density up to z=2 



translates into a stellar mass growth for the quiescent pop- 
ulations of ~ 0. 065 M© / Mpc 3 /yr, which is more than 3 
times faster than the evolution at z < 1.2. In contrast to 
z < 1.2, where most of the galaxies seem to be in place, 
at z > 1.2, galaxies are still in an active phase of mass as- 
sembly. In particular, the evolution of the quiescent popu- 
lation suggests that we are observing the epoch when , for 
the first time in the history of the universe, a large num- 
ber of active galaxies are ending their star formation ac- 
tivity and start to build up a quiescent population. While 
the mechanism acting in this process is not clear, gas ex- 
haustion, merging or other effects, this build-up happens 
a few Gyr after the peak of the cosmic SFR (Hopkins and 
Beacom, 2006). An interesting related information is the 
similar evolution followed by morphologically selected el- 
liptical galaxies (Franceschini et al., 2006; Abraham et al., 
2007). If not by chance, this coincidence could suggest that 
the build-up of the quiescent sequence is closely followed 
or preceded by a morphological transformation. 
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